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Abstract 

Models where relativistic jets from merger or accretion induced collapse 
in compact binary systems produce cosmological gamma ray bursts (GRBs) 
also predict that GRBs are accompanied by delayed emission of high energy 
photons, TeV neutrinos. X-rays, optical photons and radio waves. Their 
emission mechanisms are similar to those responsible for their emissions by 
blazars during gamma ray flares. The recently observed afterglow of GRB 
970228 in X-rays and optical photons is well accounted by such models. In 
particular, the predicted power-law spectrum, time delay which is inversely 
proportional to frequency and a power-law decay with time which is energy 
independent, are all in excellent agreement with the observations. 



Introduction 



The origin of gamma ray bursts (GRBs), which were discovered more than 
25 years ago, has been a complete mystery (e.g., Fishman and Meegan 1995). 
Their observed isotropy in the sky, deficiency of faint bursts and the lack of 
concentrations towards the Galactic center, in the Galactic disk and in the 
direction of M31 strongly suggest (e.g., Briggs 1995) that they are cosmo- 
logical in origin (Prilutski and Usov 1975; Usov and Chibisov 1975; van den 
Bergh 1983; Paczynski 1986, Goodman 1986; Goodman, Dar and Nussinov, 
1987; Eichler et al. 1989) but no conclusive evidence for that was available. 
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However, this have been changing dramaticly since February 28. On Feb. 
28.458 UT, eight hours after the gamma ray detectors aboard the Compton 
Gamma Ray Observatory, Ulysses space craft. Wind spacecraft and Bep- 
poSAX satelhte have detected the gamma ray burst GRB 970228 within a 
narrow error box (Costa et al. 1977a; Palmer et ah 1997; Hurley et al. 1997), 
the BeppoSAX satellite detected in the 0.5 — 10 keV energy range an X-ray 
counterpart, i.e., a fading X-ray source at the same position which was not 
known before (Costa et al. 1997b). The field was observed again on March 
3.734 by BeppoSAX and on March 7.028 -7.486 by ASCA in the 2-10 keV 
range (Yoshida et al. 1997) and the source was detected at the same position 
at a flux level lower by a factor ~ 20 and ~ 35, respectively. An optical 
fading remnant was also detected 21 hours after the detection of the GRB 
(Groot et al. 1997) with magnitudes my = 21.3 and mj — 20.6 in the V and 
I bands, respectively. The fading optical GRB remnant was acquired by the 
Hubble Space Telescope (HST) on March 26 which resolved it into a point- 
like source at 25.7 magnitude (Shau et al. 1997a) and an extended object (a 
faint galaxy ?). It was detected again by HST on April 7 with mv = 26.0±0.3 
and mi — 24.6 ± 0.3 (Shau et al. 1997b). In the R band, where the relative 
contribution of the faint galaxy is expected to be larger, the remnant was 
detected with a slower decline, i.e., with ~ 24 between March 6.32 -13.0, 
by Keck II, INT, Palomar 5m and ESO NTT, and with tur ^ 24.9 ± 0.3 
around April 6 by Keck II (e.g., Metzger et al. 1997). The identification of 

the site of GRB 970228 with a very faint galaxy places it at a cosmological 
distance (although no absorption features that could allow redshift determi- 
nation have been reported). Cosmological distances {D ~ lO^^cm) imply that 
GRBs have enormous energies, E ~ 10^^~^^AQ erg, during short periods of 
time, where AQ is the solid angle into which their emission is beamed. This 
enormous energy-release over a short time suggests that merger or accretion 
induced collapse (AIC) of compact objects (neutron stars and black holes) 
in close binary systems produce cosmological GRBs (Blinnikov et al. 1984; 
Goodman et al. 1987; Eichler et al. 1989). Relativistic beaming is required 
then in order to avoid self absorption through 77 — > e'^e". Relativistic fire- 
balls (Cavallo and Rees 1978, Paczynski 1986; Goodman 1986) or relativistic 
jets (e.g., Shaviv and Dar 1995a, 19996) therefore were invoked in order to 
explain GRBs. Observations and theoretical considerations favor relativistic 
jets as their drivers. Observations seem to indicate that all systems which 
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are powered by mass accretion onto a compact object release a large fraction 
of their accretion power in kinetic energy of highly relativistic jets. These 
relativistic jets are observed to produce very efficiently photons whose ener- 
gies extend from radio waves all the way up to TeV 7-rays. The required 
energy release is more moderate in the jet geometry. Therefore, Shaviv and 
Dar (1995, 1996) developed a model for cosmological GRBs where highly 
relativistic hadronic jets in dense stellar regions upscatter stellar light into 
gamma rays. 

Hadronic collisions of the relativistic jet particles with gas targets (stellar 
wind, bloated stellar atmospheres, molecular cloud) can produce high energy 
gamma rays, high energy neutrinos, and high energy electrons and positrons. 
It was suggested that these high energy electrons and positrons produce, 
through synchrotron emission and inverse Compton scattering, afterglows in 
X-rays and optical photons of 7-ray flares from blazars (Dar and Laor 1997) 
and of GRBs (Shaviv and Dar, 1996; Dar 1997). In this letter we show that 
emissions from a jet can account very well for the observed afterglow of GRB 
970228 in X-rays and optical photons. 

GRBs From Accretion Jets 

The formation, of highly relativistic jets are not well understood. However, 
the final accretion rates and the magnetic fields which are involved in mcrg- 
ers/AIC of compact stellar objects are much larger than those encountered 
in active galactic nuclei (AGN), microquasars and high mass compact X- 
ray binaries. Therefore, it is very hkely that highly relativistic jets are also 
produced in mergers/AIC of compact stellar objects. In view of the un- 
certainties in modeling mass ejection driven by merger/AIC, rather Shaviv 
and Dar (1996) assumed that NS/AIC produce jets with typical Lorentz fac- 
tors, r ~ 10^, which are required by GRBs, beaming angles Afl > tt/F^ 
similar to those observed/estimated for AGN and microquasars, and ejected 
mass AM ~ {dM/dQ)AO, < lO^^M© (because the released kinetic energy 
is bounded by Ek — FAMc^ < Mqc^). The natural birth places of close 
binary systems are dense stellar regions (DSR) such as galactic cores, star- 
burst regions and super star-clusters. These DSRs are rich in gas and in very 
luminous stars which emit very strong stellar winds (e.g., Tamura 1996). 
They have rather large photon and gas column densities, A^-y and Np, respec- 
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tively, larger than Resonance upscattering of interstellar photons 

in such DSR can explain the fluence, duration complex light curves and spec- 
tral evolution of GRBs, provided that the particles in the accretion jets are 
injected with approximately a simple power-law distribution of Lorentz fac- 
tors, dn/dV ~ with 2 < a < 3, above a minimum value Tj^ ~ 300 
(Shaviv 1996; Shaviv and Dar 1996). 

In fact, for a free ballistic expansion the ejected jet is optically thick to 
interstellar photons if AM > mA-R^AQ/cr 10~^MqAQ, where i? ~ 5 x 
lO^^cm is the typical size of DSR and a ~ lQ~^^cm^ is the typical resonance 
scattering cross section of optical photons from typical mergers/ AIC hadronic 
jets. Therefore, the ejected jet acts as a highly efficient relativistic mirror 
which upscatters stellar photons along its trajectory into beamed 7-rays. The 
typical photon fluence and energy of GRBs produced by merger/ AIC in DSR 
are 



The observed emissions of GeV 7-rays (e.g, Thompson et al. 1995) and TeV 
7-rays (e.g., Kerrick et al. 1995; Gaidos et al. 1996) from blazars are usually 
interpreted as produced by inverse Compton scattering of highly relativistic 
electrons in their jets on soft photons, internal or external to the jet (Dermer 
and Schlickeiser 1993, 1994; Blandford and Levinson 1994). However, the 
kinetic energy of relativistic jets that are made of normal (hadronic) mat- 
ter resides mainly in protons and nuclei. Therefore, the above mechanisms 
seem unable to convert enough kinetic energy of hadronic jets into high en- 
ergy photons and explain the observed large fluxes of GeV and TeV 7-rays. 
Hadronic jets, however, can convert a significant fraction of their kinetic en- 
ergy into GeV- TeV 7-rays through pp — > 7r°X; 7r° 27 in collisions with gas 
targets of high column density along the line of sight, which are present near 
accreting massive black holes (Dar and Laor 1997). Gas targets with high 
column densities (stellar winds, bloated atmospheres) can also be present 
near accreting stellar compact objects or in star-burst regions and galactic 
cores where NS merge/AIC are expected to be highly enhanced. They may 
convert a significant fraction of the jet kinetic energy into high energy 7-rays. 



^ aN^{Ve/mAC^)EK ~ 10^^^^ AQ erg. 



(1) 
(2) 



High Energy Gamma Rays 
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The cross section for inclusive production of high energy 7-rays with a small 
transverse momentum, cpr = Et < 1 GeV in pp collisions (e.g., Neuhoffer 
et al. 1971; Ferbel and Molzon 1984) is well represented by 

^ (l/27r£;o^)e-^-/^° (3) 



where E is the incident proton energy, ~ 35 mh is the pp total inelastic 
cross section at TeV energies, Eq ^ 0.16 GeV and f^{x) ~ (1 — xY / ^/x is 
a function only of the Feynman variable x = E^/E, and not of the separate 
values of the energies of the incident proton and the produced 7-ray. The 
exponential dependence on Et beams the 7-ray production into 9 < Et/ E ~ 
l/6r along the incident proton direction. When integrated over transverse 
momentum the scaled inclusive cross section becomes a~^da/dx ^ f-yix). If 
the protons in the jet have a power-law energy spectrum, dFp/dE pa AE~°', 
then, because of Feynman scaling, the produced 7-rays have the same power- 
law spectrum: 



dF^ dF„ da 



/CO rllH H(j 
^ ^-^dEp^Npa^J^AE--, (4) 

where Np is the column density of the target and = Jq x"'^^ f^(x)dx. Gas 
targets with large column densities, Np > 10^'^ cm^^, are present in/near 
galactic cores and in star burst regions. For instance, such regions seem to 
be very enriched in variable luminous blue supergiants (e.g., Tamura 1996). 
Such stars whose luminosity exceeds ~ IO^Lq, also eject mass at a rate of 
up to Lm ~ 1O~*^M0 yr~^ with an outflow speed of f ~ 500 km s^^. If 
such a star is located at a distance i?* ~ 5 x 10^^ cm from the explosion 
and at an angle 6** < l/F relative to the jet, its wind produces an effective 
gas target inside the beaming cone with an average column density of Np fa 
Lm/2mpvR^9* ~ 10^^ cm~^, which can convert ~ l^f^^^erg of the jet energy 
into high energy 7-rays. 



High Energy Neutrinos 

Hadronic production of photons in diffuse targets is also accompanied by 
neutrino emission through pp — > tt^ ^ l^^^n ; A*^ ~^ e^i/^z/g- For a proton 
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power-law spectrum, dFp/dE = AE^" with a power index of « ~ 2, one 
finds (e.g., Dar and Shaviv 1996) that the produced spectra of 7-rays and 
i/^'s satisfy 

dF^JdE ^ NpGinh^AE-'^ ^ 0.7dF^/dE. (5) 

Consequently, we predict that high energy 7-ray emission from GRBs is ac- 
companied by emission of high energy neutrinos with similar fluxes, light 
curves and energy spectra. At energies above TeV, 7 rays from distant 
GRBs are strongly attenuated by the infrared cosmic background radiation 
(e.g., Stecker et al. 1993). However, the universe is transparent to TeV 
neutrinos. Therefore, detection of TeV neutrinos can be used to confirm 
the hadronic nature of jets and to detect distant GRBs in TeV emission. 
The number of z/^ events from a GRB in an underwater /ice z/^ telescope is 
SNaT J Rn{da^^/dEfj){dFy/dE)dEndE, where S is the surface area of the 
telescope, A^^ is Avogadro's number, T is the duration of the GRB, a^^ is 
the inclusive cross section for i/^p /iX, and i?^ is the range (in gm cm"'^) 
of muons with energy i?^ in water/ice. For a GRB with F^ ~ 10~^ erg cm~'^ 
above E^ = 1 TeV and a power index a = 2, we predict 3 neutrino events in a 
1 fcm^ telescope. If the observed GeV gamma ray emission from GRB 940217 
(Hurley et al. 1994) extended into the TeV region it could have produced 
hundreds of events in a 1 km^ neutrino telescope. 

X-Rays, Optical Photons and Radio Waves 

The production chain pp ^ ti^ ^ fi"^ —>■ that follows jet-gas collisions 
enriches the jet with ultrarelativistic electrons. Due to Feynman scaling, 
their differential spectrum is proportional to the emitted high energy 7-ray 
spectrum, 

dFjdE NpainhAE-'^ ^ OAdF^/dE. (6) 

If they encounter a perpendicular magnetic field they cool radiatively by 
synchrotron emission and inverse Compton scattering and their initial energy 
decreases according to F = Fq/ (1 + t/ti/2) with ti/2(Fo) = 3mc/4aTFo(?v,B + 
Uy), where ax is the Thomson cross section and Ub = B'j^/Sn and Uy are 
the energy densities of the transverse magnetic field and the radiation field, 
respectively. The acceleration, the coUimation and the little deflection of jet 
particles by interstellar magnetic fields, as well as the emission of polarized 
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radiations by rclativistic jets, indicate that strong magnetic fields are present 
within jets. Similarly, strong radiation fields may also be be present within 
jets due to the beamed emission of synchrotron radiation (inverse Compton 
scattering from synchrotron photons of energy e can upscatter them to energy 
~ (4/3)rge in the observer frame). Consequently, for a > 2, radiative 
cooling by inverse Compton scattering and synchrotron emission changes the 
initial power-law energy spectrum of the ultrarelativistic electrons into 

dFjdE ^ NpaiJeAE-'^il - t/h/2{E))''-^. (7) 

Thus, for a > 2 the electron spectrum cuts off at ~ SrUeC^ /4:{uB + Uj)aTct 
and the photon emission cuts off at a corresponding frequency. For t <S ti/2 
the electron spectrum maintains its original spectral form ~ E~°'. Syn- 
chrotron emission and inverse Compton scattering from electrons with such 
a power-law spectrum produce a power-law spectrum: 

dFjdE ~ (8) 

The radiative cooling time of jets must be very long (jets from microquasars 
reach hundreds light years, while AGN jets reach hundred thousands light 
years before disruption). However, the magnetic (and the radiation) energy 
density within the ejected jets that expand ballisticly decreases like ub — 
TT ~ i? ~ t ^ where R is the distance from the ejection point and t is 
the time in the observer frame. Maximum emission of synchrotron radiation 
by electrons with energy Ee occurs at photon energy E^ ~ B±E'^. Photons 
with energy E^ at time t are produced mainly by electrons with energy Ee ~ 
{E^/B±y/'^ ~ E^t^/"^. Because the cooling time of jet particles is relatively 
long, the electron spectrum is proportional to E~°' and consequently the 
decay of the afterglow is independent of frequency and is given approximately 

by 

d^FjdEdt ~ £;-("+l)/2i-3{a-l)/4_ (g) 

Eq. (9) also follows from the fireball model of Wijers et al. 1997. It is 
valid until inverse Compton scattering of the microwave background photons, 
which produces the same power-law spectrum but decays very slowly with 
time, takes over. The model also predicts a time lag td in the appearance 
of the afterglow in different energy bands, which is inversely proportional 
to photon energy. This time lag is due to the difference in times of flight 
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of electrons with different energies along the jet, as seen by the observer, 
which is proportional to l/Fg, i.e., td ~ The measured afterglow in 

X-rays and optical photons of GRB 970228 at different observations times 
during the first 38 days are very well described by eq.(9) with a ~ 2.6, i.e., 
by Ed?F^/dEdt ~ E-^H-^-'^ (see, e.g., Wijers et al. 1997). 

Conclusions 

If GRBs are produced by merger/AlC jets in star burst regions or galactic 
cores of distant galaxies, they are also accompanied by delayed emission (af- 
terglow) of high energy photons, TeV neutrinos. X-rays, optical photons and 
radio waves. Their predicted spectral and temporal behavior are correlated 
and follow simple power-laws, as described by eq.(4) for high energy gamma 
rays, by eq.(5) for TeV neutrinos, and by eq.(9) for X-rays, optical photons 
and radio waves, respectively. The predictions are in excellent agreement 
with the observed afterglow of GRB 970228 in X-rays and optical photons. 
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